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Abstract 



A portion of 20 minutes of the GPS signals collected during the Bridge 2 experimental 
campaign, performed by ESA, have been processed. 

An innovative algorithm called Parfait , developed by Starlab and implemented within 
Starlab's GNSS-R Software package STARLIGHT (STARLab Interferometric Gnss Toolkit), 
has been successfully used with this set of data. 

A comparison with tide values independently collected and with differential GPS processed 
data has been performed. We report a successful PARIS phase altimetric measure of the 
Zeeland Brug over the sea surface with a rapidly changing tide, with a precision better than 
2 cm. 



1. The GPS reflected fleld retrieved by GPS processing 



1.1. Introduction 

As reported in |2^, during the last decade many experimental camp aigns to recover GPS 
reflected data have been successfully organised (a partial list in Table (|l.l.|) ). 

The data which have been used in the work described in this paper were collected during 
the Bridge-2 experiment |po|| . 

This article covers the following issue: 

• Analysis of the direct and reflected signals and EM field extraction. 

• Analysis of the reflected field behaviour. 

• PARIS phase altimetry. 

• Analysis of the altimctric results. 

1.2. Phase retrieval 

The importance of retrieving the phase of the EM field backscattered by the sea surface, 
in an Earth observation application, is related to the possibility of accurately estimate the 
delay of the reflected signal w.r.t the direct one, i.e. for altimctric purposes as well as (p^) 
is related to the sea surface conditions. 

Clearly, in order to collect the phase of the EM fleld, the complex signal has to be considered. 
In the current experiment, only the in-phase component of the signal is sampled and stored. 
It is necessary than to generate locally a quadrature component. This is done as presented 
in Figure (|^). Let us represent the signal received at the antenna input ^ as 

s{t) = C{t)D{t)cos{{ujL^ + uJd)t) (1) 

where C{t) represent the C/A code, D{t) the navigation code, luli the LI carrier frequency, 
and ujd the actual Doppler frequency. After modulation with a local oscillator of frequency 
ujli —ojif and low-pass filtering, the signal will have a residual carrier at tUd+uiiF- This signal 
is then mixed with to phasors at frequency lujf + i^d, where ujd is the available estimation of 
the Doppler frequency for the satellite under investigation. These two phasors are relatively 
delayed by 7r/2 radians, i.e. they are in quadrature. The result of such a mixing is a signal 
with two frequency component: Aw^, the error in the Doppler frequency estimation, and 
2(lojf + ijJd) — AiUd, a spurious frequency which is not possible to filter out and that is 
considered as a noise component. 

With the assumption that the navigation bit is constant during the integration time 
(always true if = 1 ms, in the other cases it is true if the code is aligned with the navigation 
bit and the coherent integration time is a sub-multiple of 20 ms) and that the component 
of the signal at higher frequency does not correlate with the clean replica of the C/A code 
(thanks to the PRN code properties) , the i-th sample of the correlation between the in-phase 
component of the signal and the clear replic can be written as 

^ iNe 

I, = -DkR, cosiAiVdJk) (2) 

{i-l)NE + l 

and for the in-quadrature component: 

Q, = -DkR, J2 siniAujdJk) (3) 

{i~1)Ne+1 

*This represent only the C/A code part. The P code part can bo neglected thanks to the subsequent 
correlation of the signal with replicas of the C/A code. The two code are in fact orthogonal. 



Main 
Author 


Affihati 


3rDatc 


-Note 


Garrison M 


JPL 


1996 


A normal receiver was used. Demonstrated 
reception/tracking of reflected signals. Con- 
cluded that more complex receiver is needed. 


Martin-Neira 01] 
Caparrini 


ESA 


1997 


PARIS altimetric concept applied GNSS-R 
collected from a bridge. C/A code only could 
be used for correlation leading to an altimetric 
accuracy in the order of 3 meters (1% of the 
chip length). 


Garrison 


JPL 


1997 


Widening of correlation function in rough seas 
demonstrated. 


Komjathy |15| 




1998 


Aircraft experiments. 3-5 km altitude. 


LaBrecque 




1998 


The first spaceborne observation of GPS sig- 
nals reflected from the ocean surface. 


CardcUach g@ 

Ruffini 

Garrison 


lEEC 


1999 


Balloon experiment. Successful detection of 
reflected signals at 38 km of height. 


CardcUach M 
Rufhni 


lEEC 


1999 


Aircraft experiment. 


Ruffini 1) 
Caparrini 




2000 


The set of GPS-R data collected from a 
drilling platform has been analysed. 


Armatys [0 




2000 


Wind speed and directions obtained from re- 
flected GPS signals are compared to the Sea- 
Winds scatteromcter on-board QuikSCAT. 


Garrison 


JPL 


2000 


With GPS-R airborne data, retrieval of the 
wind speed with a bias of less than 0.1 m/s 
and with a standard deviation of 1.3 m/s. 


Zavorotny |^ 




2000 


Comparison of experimental and theoretical 
waveforms. 


Zuffada H 


JPL 


2000 


Lakeside experiment, with an almost local pla- 
nar surface (no roughness). 


Martin-Neira [Ij |lj 


ESA 


2000 


The experiment was designed to test some key 
issues in the PARIS Interferometric Proces- 
sor (PIP) concept. The PIP concept is based 
on the use of dual-frequency carrier measure- 
ments to exploit the correlations in the scat- 
tered signals at similar frequencies. 


Martin-Neira |^ 


ESA 


2001 


The experimental campaign which is the ob- 
ject of this work. 


CardcUach y y]j5|[y 

Germain 

Caparrini 


lEEC 
STAR- 
LAB 


2001 


GPS-R data collection from airborne plat- 
form. Campaign performed within the 
ESA/ESTEC project "PARIS Alpha". Data 
processed under ESA/ESTEC projects 
"PARIS Alpha" and "OPPSCAT 2". 


STARLAB 


STAR- 
LAB 


2002 


GPS-R data collection from airborne plat- 
form. Data processing currently performed 
within the ESA/ESTEC project "PARIS 
Gamma" . 



Table 1: Representative GNSS-R experiments (1996-2002). 



Considering that during an integration time interval the value of Ao^d^ is constant, joining 
the two signal components and performing the summation, the complex i-th sample of the 
correlation writes 



1 „ sin I^^Te] 



C, - -DkR^e-^^'^'^^^^'-^^e-^'^'^'^^ >- f (4) 



sm -^T. 



Our final goal is to track the value of the carrier phase. To this end, the delta-phases 
obtainable from equation (jj) are accomulated as follows: 

- 0, = $,+1 = (log^) = -^uJdTE. (5) 

This equation is true if Awc;.^^ « Awtj. , which is a good approximation since the time dur- 
ing which this variation is measured is the coherent time of integration, i.e. in the order of 
milliseconds. 



The main advantage of using this delta-phase approach to phase tracking is that we have 
to deal with differential value, that is smaller values than the full phase values. This allows 
an easy detection of the tt radiants change of phase due to the presence of the navigation bit. 

Figures from ||to (Q), used to explain these concepts, refers to the processing of another 
set of GPSR data (the Casablanca Experiment). In Figure (||) the histogram of the delta- 
phases is shown. The x-axis represent cycles and the y-axis just arbitrary units. The most 
part of (5-phase values is clearly concentrated around zero. Another considerable amount of 
values is allocated just before ±7r. These values represent in fact small values to which ±7r 
radiants have been added, for the occurrence of a navigation bit. 

In Figure (H) both the phase with and without the navigation bit effects are plotted. In 
order to derive geophysical parameters from GNSS-R, the phase has to be considered without 
any Doppler contribution. In Figure (|^) this phase for the direct and the reflected signal are 
shown. The effect of the reflection on the sea surface is clearly visible in the huge variations 
present in Figure ( 4(b) ) with respect to Figure ( 4(a) ). In Figure (||) and (^) the amplitude 



magnitude and respectively the complex vector of the direct and reflected fields are shown. 
In Figure (j^) a simulation of the LI GPS field after reflection is shown. The parameter of 
the simulation are chosen according to the Casablanca experiment characteristics. 



2. A PARFAIT approach 
2.1. Introduction 

It is know that the major inconvenient while tracking the reflected field phase is the occur- 
rence of fadings. Even though in the Bridge 2 experiment, due to the relative smoothness of 
the sea surface [^^, these problematic events arc not so frequent, nonetheless, since a single 
event can completely destroy the phase information, some countermeasures must be taken. 
For example (see [pO|), it is possible to inject in the system (during the fading events) a 
model-based phase. Another possibility is not to track the reflected field phase at all but to 
track the direct field minus the reflected fleld phase, which is really what PARIS concept is 
about. 



Let us see more in detail why we should not try to track the phase at all. First, a 
distinction should be made. Using the accumulated phase for altimetry presupposes two 
difficult tasks: 

• Tracking the phase (hard due to bandwidth and fadings). 



• (If tracking is possible) use the tracked phase for altimetry. 



Note that the second point is often overlooked. Yet we know that the accumulated phase 
acts like a random walk. The average accumulated phase is very badly behaved (essentially 
meaningless). 

Tracking the phase is important if the signal of interest rotates fast as compared to noise- 
induced rotation. This is definitely not the case in the bridge phase or in any other coastal 
application of the technique. In fact, it is never true for PARIS applications over the oceans. 
It is also very difficult to do: errors accumulate fast, as above mentioned for example, during 
fading events. 

Even more important, the accumulated phase is meaningless. As we have shown in previ- 
ous work, the phase wanders around, jumping to different winding number kingdoms. This 
would appear to be a disaster, but it is not: what we would like to get, if the geophysical 
signals of interest move slowly enough, is the average field. This exists even if the phase 
wanders and wanders. 

In the bridge experiment the sea state was very smooth. In such a case, tracking the 
phase may still be useful but in general, this is definitely not true. 

2.2. Using the direct minus reflected combination 

What we are after is the difference of direct and reflected signals. Note that the same lo^ 
will be used for direct and reflected — this simplifies things. We can then return to equation 
(|5|), and write 

Im [In - In Cf )] - d^f = - At^^-^T^j. (6) 

Note that we are measuring directly the interferometric field. Let us now address some of 
important points. 

2.3. Why we don't need to worry about emission versus reception 
time diflFerences 

The difference between direct and reflected paths is much less than a microsecond (300 light- 
meters). Even for a fast moving GPS satellite, this is essentially zero. 

This is the essential measurement we arc making on reception of a GPS signal: 

cM = \\xB{tR) - XE{tE)\\ ~ \\xR{tR) - XEitn) - VE5t\l (7) 

where 6t is the time it take from emission to reception. This is of the order 2x10^ light 
meters (about 0.07 seconds). Definitely significant: the emitter will move about 300 me- 
ters in this time. What does it mean? Well, since we will work with the direct minus 
reflected travel times, the effect is null. Recall that from a ground platform we essentially 
have ^Db. ~ 2/isina, where a is the elevation. In 0.07 seconds, the change of elevation is 
negligible. 

On the base of these considerations, a new approach to PARIS altimetry is described in 
the following sections. It is called PARis Filtered-field Altlmetric Tracking technique (PAR- 
FAIT technique) and it proved to be a very robust and precise processing method. 



3. The PARFAIT processing 

In this section the basis of the parfait processing arc presented. The parfait pro- 
cessing is an innovative processing for GNSS-R data developed by Starlab, as part of the 
STARLIGHT GNSS-R software package. In the frame of this CCN, the application of such 
a processing will lead to the estimation of the height of the antennas of the Bridge 2 Experi- 
ment over the sea surface. 



Simple geometrical considerations lead to the following equation which relates the height 
of the receiver over the reflecting surface (considering the same height for the upward looking 
antenna and for the downward looking antenna) with the delay measured between the direct 
and reflected GPS signals: 

delay PRN {t) ■ c = 2h (t) sin {epRN (*)) + doff set (8) 

where delayppiq{t) is the measured delay between the two signals at time t, hit) is the height 
of the bridge at time <, epB,N{t) is the elevation at time t of the GPS satellite with a specific 
PRN number, and doff set is the hardware-induced delay, considered to be constant in time 
(c is the light speed). From this equation it is evident that a first estimation of the height 
of the receiver can be easily done with a linear fit of the measured delay with respect to the 
sinus of the elevation angles of each satellite. The delay between direct and refiected signal, 
though, can be measured through two quantities: C/A code and carrier phase. In the case 
of carrier phase, equation (||) must be rewritten as follows 

delaypRN (t) ^ 2h (t) sin {eppN (t)) + Nprn>^ + doff set (9) 

where this time delay ppN {t) is already in meters and A is the carrier wavelength. In other 
words, the equation of each satellite contains an unknown parameter NppN. In order to use 
all the satellites for one height estimation, it is then necessary to estimate also NppN, i-e. to 
solve the ambiguity problem. In practise, the following procedure was used. All the measured 
delay ppN(t) have been put in an unique vector, with an additional additive parameter bpuN 

{delayi{ti), delayi{t2), ((ieZai/2(ii) + ^2), ■■• 

{delay2{t2) + b2) , {delayp{t„i) + bp)) . 

The corresponding vector for the sine of the elevation angles was then built. Using these 
two vectors, a linear fit has been performed, resulting in a certain value for the norm of the 
residuals of the fit. Finally, considering this norm as a function of the parameters bppN, the 
optimum value bpp^ that minimises the norm has been found. For this bppN one obtains 
the best interpolation, i.e. the best value for a first estimation of the height of the bridge. 
This estimation of the height is nonetheless quite rough. Let us see why and how to improve 
this estimation. First of all, remember that the term on left-hand side of equation (^) can be 
written as (see also equation (^)) 

delaypRN (t) = A^f"^ (10) 

where d(j)^~^ is the phase of the (refiected minus direct) interferometric field. This field 
is generally corrupted by fadings and, during fadings, the phase of this interpherometric 
signal is completely impossible to track (basically, during these events, the reflected signal 
is not present at all). Moreover, the optimisation previously described is performed in the 
TZ domain, without correctly approaching the ambiguity problem which is inherently to be 
solved in the integer domain. 

To face the fading problem, we propose to filter the interferometric field. This filtering should 
be long enough to eliminate fading problems and short enough to let the signal of interest 
pass through. Consider the usual equation 

Xd(t>f"°{t) = 2h{t) sin{epRN{t)) + doff set. (11) 

The value of dtpf'^^ should not change for more than a fraction of 2tt in the time duration of 
the filter. This maximal time depends clearly on the elevation angle of the satellite and, just 
slightly, form the tide motion. In the case under examination the maximum filtering time is 
around 10 seconds. In other words, in 10 seconds, at least for one satellite, d4>f~^ changes of 
^ radiants. With this filter length, the fadings are not completely eliminated, even though a 
realistic estimation of the bridge height (and doff set) is now possible. To definitely override 
the fading problems, a longer integration should be performed, without killing the signal 



{d(f>^ ~ ) we wish to measure. To this end, it is possible to counter-rotate the interferometric 
field, using the first guess of the bridge height as explained in the following. Consider the 
interferometric field as 

E.ntevS oc exp{d^f-°} = exp{2/i(t)... (12) 

sin(epflAr(t)) + (io//set} (13) 

and assume a first guess ft.;, for the bridge height (and doff set for the offset). Than the 
interferometric field is down- converted or counter-rotated using the bridge height initial guess 

EZ?'erf « exp{2 h it) sin (ep^ (t)) + ... (14) 

doff set} ■ exp -2 hb (t) sin (epRN (t)) + doff set (15) 

obtaining 

Etn?erf « exp{2 Shb (t) sin (epRN (t)) + Sdoffset}- (16) 

Clearly, the phase of the field in equation ( [l^ ) varies much slower than the phase of the orig- 
inal interferometric field and this allows a longer filtering time, that is, virtually, a complete 
elimination of fadings problems. 

This technique provides a fundamental cornerstone for PARIS processing from air and space- 
borne applications. 

The equation which relates the phase delay between direct and reflected signal, the satellite 
elevations and the 6hi, (i.e. the error between the first guess of the bridge height and the real 
value) turns out to be 

Ad0f-^ (t) = 2 Shb (t) sin {epjiN (i)) + NprnX + Sdoffset (17) 

This is the new equation to be used to fit the straight line and infer the height of the bridge, 
with respect to the first guess used to counter-rotate the field. As already said, the equation 
relating the interferometric phase and the bridge height is known up to a multiple of wave- 
length (equation (^)) or, in other words, there is an ambiguity due to the periodic behaviour 
of the phase. In order to correctly solve this ambiguity problem, a search is performed in the 
space of the n-uples of integers and the one that produce the linear fit with smallest residue 
is assumed as the true one. It is important to point out that the space of the n-uples to be 
spanned in order to remove the ambiguity is drastically reduced by the filtering of the field 
previously described. In fact, since the only residual movement in the filtered data is the 
tide, the single filtered phase history, will all lie in a narrow height interval. This interval is 
determined basically by the correctness of first guess for the bridge height. In other words, 
if the guess is within ± half meter, the n-uples subspace to be scan can be limited to those 
n-uples whose components belong to the interval [—3,3], centred on the first guess of the 
n-uple, obtained from the real solution. 

Another possibility to reduce the cardinality of the subspace of the n-uples to check, is to 
consider that satellites with similar elevation angles cannot have very different integer ambi- 
guities. 



4. First altimetric results 

The PARFAIT analysis described in the previous section, has been applied to the first 10 
minutes of the Bridge2 data, part Al and to the first 10 minutes of the part A2. The following 
step have been performed accordingly: 

• the EM fields, direct and reflected, have been computed through the GPS correlation 
process; 



PRN 


elevation 


mean SIMK 
(direct) 


mean SIV K 
(reflected) 


14 


17° 


29.4 dB 


25.0 dB 


25 


17° 


32.0 dB 


25.8 dB 


1 


30° 


31.2 dB 


24.6 dB 


7 


38° 


33.2 dB 


29.4 dB 


11 


62° 


34.0 dB 


29.4 dB 


20 


78° 


30.4 dB 


26.6 dB 



Table 2: Visible satellites, their elevation in degrees, the 10 ms coherent integration mean SNR^g^^ 
direct and the reflected signal. 



for the 



• the field has been comrter-rotated (equation (^); 

• the field has been filtered; 

• the phase of the field has been evaluated; 

• a straight line has been interpolated to the phase histories (one for each visible satellite) 
against the elevation angle (equation (|T7|)); 

• the real n-uple of values of XNpn^ that minimise the residue of the fit is use to define 
a first guess for the integer n-uple of values of Npji^^ ; 

• the linear fit has been again performed for every integer n-uple of Npfij\j to find the 
one that minimises the residuals. 

The filter used for smoothing the field is a fiat zero-phase filter, with a length of 30 
seconds. 

This filtering has been performed for almost]^ a ll visi ble satellites (see table (^) and Figure 
The phase histories are shown in Figure ( |lO(a)| ). A straight line has be en fit through 
these phase histories, against the sinus of the satellite elevation angle (figure ( 10(b)| )). 

This fitted line gives an estimation of the bridge height of 18.61 m, an hardware bias of 
-0.81 and, as first guess for the n-uple that solves the ambiguity problem, the values [001123]. 
Now, a search in a subset of is done to minimise the residuals of the fit in the space of 
the n-uplcs of integers. The subspace considered is the one spanned by all the combination 
of integers between ±3 around the first guess. The result is the n-uple [002245] which gives 
a bridge height estimation of 18.82 m and an instrumental bias of —0.45. 

This procedure has been applied to the first 10 minutes of the part Al and of part A2 
of the data. The results are reported in Table (||) and in Figure (^l|) for part Al and in 
Table (||) and in Figure (13) for part A2. The final value of the estimated bridge height is 
considered to be the interpolated straight line through the available point after removing the 
supposed estimation bias . The use of a straight line for the interpolation is justified by the 
short period of time considered, relatively to the tide period. The standard deviation of the 
interpolated estimation w.r.t. the measured tide is of 0.35 cm for part Al and of 0.84 cm for 
part A2. The results for both periods arc plotted in figure (0). 



Fitting both parts to the tide curve, i.e. choosing the bias that minimises the standard 
deviation of the data to tide "ground truth", leads to a bias of 41.10 cm and a standard 
deviation of 1.81 cm. This bias could be due either to an error in the determination of the 
absolute value of the height of the bridge performed with the differential GPS processing or, 
partially, to some anomalies in the fiowing of the water in the vicinity of the bridge structures. 

Moreover, assuming also that the tide measured below the bridge can have a time delay 
with respect to the place were the tide is measured, the best fit (over both bias and delay) 

t Satellites outside the Zeeland Mask , as defined in |19| are discharged (see also caption of Figure Xt). 



is obtained with a delay of 3 minutes and 12 seconds with respect to the time of the tide 
data collection and with a bias of 39.13 cm. The standard deviation of the fitted data with 
respect to the tide curve is in this case of 0.893 cm. 



time 
(minutes 
from start) 


instrumental 
bias [cm] 


bridge height 
estimation [m] 


assumed 
ground truth 
[m] 


difference 
[cm] 


1 


-0.45 


18.83 


18.44 


39.06 


2 


-0.45 


18.82 


18.42 


39.13 


3 


-0.46 


18.81 


18.41 


39.78 


4 


-0.45 


18.79 


18.40 


39.80 


5 


-0.26 


18.78 


18.38 


39.29 



Table 3: Results of the bridge height estimation during the first 10 minutes of the part Al data. 

Another interesting check to be done to asses the validity of the processing is to compare 
the change in the height of the bridge between the first 10 minutes of part Al and the first 10 
minutes of part A2 as retrieved by the processing and as provided by the tide measurements. 
This comparison is shown in figure (|^). The result is absolutely satisfying: the estimation is 
in accordance with the measures within about 4 cm. 

To summarise, the proposed approach to PARIS altimetry, the parfait technique, leads 
to possibly biased but very precise estimation of the tide, 

• without the need to insert any kind of model for the phase of the reflected signal during 
fadings; 

• without rejecting too many visible satellites because of their poor SNR and/or frequent 
fadings; 



time 
(minutes 
from start) 


instrumental 
bias [cm] 


bridge height 
estimation [m] 


assumed 
ground truth 
[ml 


difference 

[cm] 


1 


-0.27 


17.54 


17.11 


42.6 


2 


-0.28 


17.52 


17.08 


44.2 


3 


-0.26 


17.47 


17.05 


42.3 


4 


-0.08 


17.44 


17.02 


42.1 


5 


-0.08 


17.41 


16.98 


43.1 



Table 4: Results of the bridge height estimation during the first 10 minutes of the part A2 data. 
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Figure 1: Starting from the in-phase component of the sampled signal, with a carrier frequency equal to the 
sum of the IF of the receiver and the Doppler frequency, two downconvertions are performed, with two phasor 
relatively delayed of it/2. After a low-pass filtering, the two obtained signals bring information about the 
amplitude of the backscattered EM field for both the in-phase and in-quadrature components. 
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Figure 2: The histogram of the carrier phase variation, measured on an integration time interval (in this 
case 1 ms). The x-axis represents cycles, while on the y-axis there are arbitrary units. It is evident the 
accumulation of (5-phase values around zero, as well as in the vicinity of ± half a cycle. 
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Figure 3: The carrier phase obtained accumulating the (5-phase according to equation {^. The stepped 
plot represent this accumulated phase as it is , i.e. without compensating for the navigation bit half-cycle 
variation. The lower curve represent the same phase after removal of this effect. The unit of the x-axis is 
milliseconds, on the y-axis is cycles. 
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(a) Direct field. 



(b) Reflected fleld. 



Figure 4: Example of tracked phase, without the Doppler contribution. The units arc milliseconds on the 
X-axis and cycles on the y axis. The integration time is 20 ms. 
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(a) Direct field. 



(b) Reflected flcld. 



Figure 5: Example of fleld amplitude time series. The units are milliseconds on the x-axis and correlation 
units in dB on the y axis. The integration time is 20 ms. 
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(a) Direct field. 



(b) Reflected flcld. 



Figure 6: Example of complex field time series. The units are correlation units on both axis. The integration 
time is 20 ms. 



Figure 7: The EM field at LI frequency, after reflection on the sea surface. This is a simulated field. The 
simulation has been performed with a software |2l[ developed by Starlab. The simulation considers a wind 
speed Uio = 3 m/s. 
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Figure 8: The blue curve is the phase of the interferometric field for PRN number 7, from minute 2 to 
minute 4 of the part Al of the Bridge 2 experiment. It is clearly visible the occurrence of a long fading (on 
the left-hand side of the plot) and of isolated phase slips. These phenomena disappear in the phase of the 
filtered interferometric field (red line). 




Figure 9: Each coloured arc represents the position of a GPS satellite from the start of the part Al of the 
experiment to the beginning of part A2 plus 10 minutes. The PRN number is closed to the beginning of 
the data. The green mask represent the area where the GPS signal reflections are supposed to be free of 
shadowing phenomena due to the bridge structure, and therefore only the satellite within this mask can be 
taken into consideration for PARIS processing. The bold parts of the lines represent the first and the second 
10 minutes periods. 
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(a) In this figure, the reflcctcd-minus-dircct phase delay 
for each PRN is plot with Np^]\j = 0. 
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(b) In this figure, the refiected-minus-direct phase delay 
for each PRN is plot with = {001123}. 



Figure 10: Each coloured spot represents the reficcted-minus-direct phase delay VS satellite-elevation for a 
certain PRN number. 
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(a) The upper line is the bridge height estimated with the 
procedure explained in this chapter. The bottom line is 
the bridge height according to the available tide measure- 
ments and the GMV calibration of the absolute height 
4:57:27 UTC, using GPS code differential processing 
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(b) In this plot, the estimations of the height (red *) and 
the actual height (line) are shown, after adding to the 
actual height the mean of the values of the last column 
of Table (|). 




Figure 11: The bridge height estimation during the first 10 minutes of the part Al of the data. 
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Figure 12: This line represents the difference of two differences: the difference of the measured bridge heights 
during the collection of the first part of data and those measured during the second part minus the same 
difference but calculated with the estimated data. The result confirms again the effectiveness of the PARFAIT 
processing. 
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(a) The upper line is the bridge height estimated with the 
procedure explained in this chapter. The bottom line is 
the bridge height according to the available tide measure- 
ments and the GMV calibration of the absolute height 
at 14:57:27 UTC, using GPS code differential processing 
0. 
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(b) In this plot, the estimations of the height (red *) and 
the actual height (line) are shown, after adding to the 
actual height the mean of the values of the last column 
of Table (§). 



Figure 13: The bridge height estimation during the first 10 minutes of the part A2 of the data. 
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Figure 14: The solid line is the distance between the up-looking antenna and the sea surface, according to 
the available tide measurements and the GMV calibration of the absolute height at 14:57:27 UTC, using GPS 
code differential processing [pO|. The green dots are the estimated values, after removing the bias. 



